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Abstract 

Current matrix elements and observables for electro- and photo-excitation 
of baryons from the nucleon are studied in a light-front framework. Relativistic 
effects are estimated by comparison to a nonrelativistic model, where we use 
simple basis states to represent the baryon wavef unctions. Sizeable relativistic 
effects are found for certain transitions, for example, to radial excitations such 
as that conventionally used to describe to the Roper resonance. A systematic 
study shows that the violation of rotational covariance of the baryon transition 
matrix elements stemming from the use of one-body currents is generally small. 
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1 Introduction 

Much of what we know about excited baryon states has grown out of simple nonrel- 
ativistic quark models of their structure. These models were originally proposed to 
explain the systematics in the photocouplings of these states, which are extracted by 
partial-wave analysis of single-pion photoproduction experiments. This method can 
only give us information about baryon states which have already been produced in 
-nN elastic scattering, since a knowledge of the coupling constant for the outgoing 
Ntt channel is necessary for the photocoupling extraction. Photoproduction experi- 
ments have also tended to have limited statistics relative to the ttN elastic scattering 
experiments. 

The traditional theoretical approach is to describe the nucleon and its excitations 
using wavefunctions from nonrelativistic potential models, which describe baryons as 
being made up of 'constituent' quarks moving in a confining potential. This potential 
is provided by the interquark 'glue', which is taken to be in its adiabatic ground state. 
The quarks interact at short distance via one-gluon exchange. The electromagnetic 
current is also calculated using a simple nonrelativistic expansion of the single-quark 
transition operator, i.e., in the nonrelativistic impulse approximation. Not surpris- 
ingly, the resulting photocouplings and electroexcitation helicity amplitudes are frame 
dependent, with the problem becoming more severe when the photon transfers more 
momentum; current continuity is also violated. This is partly due to the nonrelativis- 
tic treatment, and partly due to the lack of two and three-body currents which must 
be present in this strongly-bound system. 

Much more can be learned about these states from exclusive electroproduction 
experiments. The photocouplings are the values of transition form factors at the real- 
photon point Q^ = (here Q^ = —q"^, where g^ is the four- momentum transferred 
from the electron); electroproduction experiments measure the Q^ dependence of these 
form factors, and so simultaneously probe the spatial structure of the excited states 
and the initial nucleons. Both photoproduction and electroproduction experiments 
can be extended to examine final states other than Nir, in order to find 'missing' 
states which are expected in symmetric quark models of baryons but which do not 
couple strongly to the Nir channel fl], 0. Such experiments are currently being carried 
out at lower energies at MIT/Bates and Mainz. Many experiments to examine these 
processes up to higher energies and Q^ values will take place at CEBAF. 

The success of nonrelativistic quark models in describing the systematics of baryon 
photocouplings does not extend to the electroproduction amplitudes. The best mea- 
sured of these amplitudes are those for elastic electron-nucleon scattering, the nucleon 
form factors. In a simple nonrelativistic model the charge radius of the proton is too 
small by a factor of almost two, and the form factors fall off too rapidly at larger 
Q^ values. It has been suggested that the problem with the charge radius is due to 
the neglect of relativistic effects 0. Similarly, the poor behavior at large Q^ can 
be improved by going beyond the nonrelativistic approximation 0, |, |, 0] and by 



enlarging the limited (Gaussian) basis in which the wavefunctions are expanded PJ. 
Similar problems exist in the description of the Q^ dependence of the transition form 



factors for the A| (1232) resonance, which are also quite well measured. Although 
the experimental information about the transition form factors for higher mass res- 
onances is limited , there are serious discrepancies between the predictions of the 
nonrelativistic model and the extracted amplitudes here also |1^, |ll], |12|. 

It is clear that, once the momentum transfer becomes greater than the mass of 
the constituent quarks, a relativistic treatment of the electromagnetic excitation is 
necessary. A consistent approach involves two main parts. First, the three-body 
relativistic bound-state problem is solved for the wavefunctions of baryons with the 
assumption of three interacting constituent quarks. Then these wavefunctions are 
used to calculate the matrix elements of one-, two- and three-body electromagnetic 
current operators. Results of model studies of the relativistic three-nucleon problem 
are available |T3|, and the problem of three constituent quarks in baryons is presently 



under investigation |I4[ (see also the comments in the final section). A useful first 
step is to consider transitions between substates in a general basis in which such 
wavefunctions may be expanded, such as a harmonic-oscillator basis. One can then 
examine successive approximations to the baryon wavefunctions with these un-mixed 
oscillator basis states as a starting point. 

There are many competing approaches to the calculation of electromagnetic cur- 
rents in relativistic bound states; we will use light-front relativistic Hamiltonian dy- 
namics |]15| here. The basis of this approach is to treat the constituents as particles 
rather than fields. It allows for an exact relativistic treatment of boosts, by remov- 
ing the dependence on the interacting mass operator of the boost generators. As a 
consequence certain rotations become interaction dependent. Within this framework, 
one can set up a consistent impulse approximation, which minimizes the effects of 



many-body currents [T^. Current continuity can be enforced, though, as we shall 



see, not uniquely, and one can explicitly evaluate the degree to which the currents 



have the required rotational covariance properties |16|. 

Other authors have applied similar models to the study of these problems; Chung 
and Coester 0], Schlumpf 0, Weber 0, and Aznaurian have studied nucleon elas- 
tic form factors in light-cone models with simple nucleon wavefunctions, and these 
models have been extended to the electroproduction amplitudes for the A(1232) reso- 



nance by Weber |T^, in a model based on light-cone field theory, and by Aznaurian 0. 
Bienkowska, Dziembowksi and Weber |]I8[ have also applied a light-cone model to the 
study of the small electric-quadrupole El multipole in A| (1232) electroproduction. 



and Weber |T9[ has also applied his model to the 5*11(1535) resonance. 

Our goal is to examine the sensitivity of these photoproduction and electroproduc- 
tion amplitudes to relativistic effects. Within the crude approximation of identifying 
un-mixed harmonic-oscillator wavefunctions with the states that they represent in 
zeroth order in the harmonic basis, we will examine the nucleon form factors, and 
the photo- and electroproduction amplitudes for the positive-parity states A| (1232), 



the Roper resonance N^ (1440), N^ (1710), A| (1600), and several negative-parity 
resonances including N^ (1535) and A^| (1520). We will show that relativistic ef- 
fects are in general appreciable; the most striking relativistic effects we have seen 
are in the electroproduction amplitudes for radial excitations such as A^| (1440) and 

A| (1600). We also provide a systematic study of rotational covariance of the baryon 
transition matrix elements. 



2 Conventions 



Much of the light-front notation is presented in Ref. |T5[. We present here the salient 
features needed to describe the calculation of current matrix elements. 

2.1 Kinematics 

A homogeneous Lorentz transformation A describes the relation between between 
four- vectors, viz., 

pM ^ p't^ = A'^^p^ (1) 

We also make use of the corresponding SL{2, C) representation, denoted by A. A 
rotationless boost is given by 

L,(A) =exp[|(A-cr)sinh'^|A|]. (2) 

Light-front kinematic variables are expressed in terms of four-vectors which are 
decomposed as a'^ = (a~, a^, a+), where a^ = a" it fi ■ a and n ■ a_L = 0. We adopt 
here the usual convention of setting fi parallel to the z axis. It is also convenient to 
define a light-front vector 

a = (a^,a+). (3) 

The SL{2, C) representation of a light-front boost is 

Lf{k) = expi^a^ In A+) exp^A^ + tA^){a^ - ta^)]. (4) 

2.2 State Vectors 

Free-particle state vectors |p/i) are labeled by the light-front vector p and satisfy the 
mass-shell condition 

p~ = {m^ + pD/p^. (5) 

They are normalized as follows: 

(pVIp/x) = {2nf6,,,6{p' - p) = {27r)%,^6{p'^ - p^)6{p'+ - p+). (6) 



The state vectors introduced above employ light-front spin. Under a light-front boost, 
p//i _ A'^j.p'^, the state vectors undergo the unitary transformation 



/+ 



[/(A)|pM) = ^^|p'/i), 



(7) 



with no accompanying Wigner rotation. These state vectors are related to this with 
ordinary (or canonical) spin via the relation the relation 



|p/^) = Y. 



l^m[P) 
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\Pf^)cDfil [Rcfip)], 



where the Melosh rotation 20 
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(9) 



2.3 Three-Body Kinematics 



Consider three free particles with masses mi, m2 and m^. We label their respective 
light-front momenta in an arbitrary frame by pi, p2 and pa. The total light-front 
momentum is 

P = Pl + P2 + P3. (10) 

Let ki, k2 and ka be the ordinary three- momenta in a frame where the total momen- 
tum is zero: 



Ek. = o. 



(11) 



i=l 



The two sets of momenta are related as follows: 



ki± 
ki3 



Pij 
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(12) 



where Mq is the invariant mass of the three-particle system: 
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Mo 



Y^u;^X^,) = ^P+P--P 



i=l 



Li=l 



Jji 



where 



co'm(k) = y/m"^ + k^. 
The Jacobian of the transformation (pi,P2;P3) -^ (P;ki,k2) is 



5(Pl,P2,P3) 



d{PMM 



PiPtptMo 



UJrm (ki)t^„2 (k2)^m3 {^?i)P^ 



(13) 
(14) 

(15) 



3 Current Matrix Elements 



To calculate the current matrix element between initial and final baryon states, we 
first expand in sets of free-particle states: 



(M'j;PV|/+(0)|Mj;P/.) 



(2vr) 



'18 



dPi / dP2 / dp3 / dpi / dp2 / dps ^ 



(16) 



{M J ;P/i|pi/iiP2/i2P3/^3) 
x(Pl/UlP2/U2P3/i3|Mj; P/i). 

Here and henceforth, an empty summation sign denotes an implied sum over all 
repeated indices. We compute only the contributions from one-body matrix elements: 

(p;/i'lP2/i2P3/^3|/^(0)|pi/iiP2^2P3/i3) 



E(p:/^:i^^(o)ip./^.)nK>.(2^)'^(p' 



^'. 



Pj 



j=l 



;i7) 



The matrix elements for the struck quark are then written in terms of the Pauli and 
Dirac form factors for the constituent quarks 

Q 



(p'/iV'-WiP/i) = Fi,(g')VM - ^K 



-F2,{Q'). 



li 



Note that this implies that these current matrix elements depend only on Q (q is 
taken to lie along the x-axis) and not on the initial and final momentum of the struck 
quark, as in the nonrelativistic model. 

The baryon state vectors are in turn related to wave functions as follows: 



(pi/ilP2/i2P3/U3|Mj; P/i) 



5(Pl,P2,P3) 



5(P,ki,k2 

xDplJ^[R,f{ks)], 



(2/^15/^2P12Ail2) (-512/^12 0/^3 PAts) 



(19) 



The quantum numbers of the state vectors correspond to irreducible representations 
of the permutation group. The spins (si2, s) can have the values (0, |), (1, |) and 
(1, |), corresponding to quark-spin wavefunctions with mixed symmetry (x^ and x^) 
and total symmetry {x^), respectively |21]. The momenta 



K, 



1 



\/6 



;ki 



2k, 



(20) 



preserve the appropriate symmetries under various exchanges of ki, k2 and ka. Note 
that kp and Ka as defined above also correspond to the nonrelativistic three-body 
Jacobi momenta. This is a definition of convenience, rather than a nonrelativistic 
approximation. The three momenta kj are defined with relativistic kinematics, and 
the use of kp and K^ accounts for the fact that only two of the kj are independent. 
It is also convenient for keeping track of the exchange symmetry of the three quarks. 
This definition allows use of the usual three-quark harmonic-oscillator wave functions 
as a basis, but there is nothing nonrelativistic about this choice. 

4 Multipole Invariants 

From the current operator /^(a;), we define an auxiliary operator /^^(x), which has 
explicit components 

I',{x) = I^x); 

/°(x) = I%x). (21) 



With these definitions, it was shown in Ref. [^ that the current matrix element 
between states with canonical spins and instant-form three-momenta could be written 
as follows: 

,(M7;Py|PfO)|Mj;P/i), 



M' 



'::^^i!^'-^>'7i^<'"'''i^'"><'^'"i*> 



X {lfiiSIJs\JfJ'j){JlJ'JlJ'j\j'lJ'') 

p' p' p' 



><yL*iPo)[L-\—)]auj7j)]lMiMUj7j),Po)]{M'j' II irsAq') \\ Hm 



where 

Po = L-\^)P. (23) 

The Wigner rotation is 

mu^ip.)^l:^(^^)U^)u^) (24) 

All dynamical information is contained in the reduced matrix element (M'j' \\ 
Iisj{q^) II Mj). The number of independent matrix elements is limited by the number 
of combinations of /, s and J which can couple together. Parity considerations provide 
the additional restriction: 

n'n(-i)'+^" = +1. (25) 



where 11' and 11 are the intrinsic parities of the final and initial states, respectively. 
Time reversal provides another constraint for the case of elastic scattering. Finally, 
the continuity equation: 

g^(M'Py|F(0)|MP^) = (26) 

further restricts the number of independent matrix elements. 

The matrix element of /^(O) between light-front state vectors is 

(My;PV|/ +(0)|Mj ;P/i) 



V P'+ V P 

X {lfiiSjis\Jfij){jfiJf^j\j'fJ^') 
p/ p/ 



^E [-(l«i/^|io) + (^"i/^|oo)] mvlsf^s) 



xFJ;(Po)[L7i(— )]c«[L^(— )]^Z^^;[P,;(Po)](M'J' H IrsAg') II Mj)(27) 



A knowledge of the reduced matrix elements {M'j' \\ Iisjiq'^) \\ Mj) is sufficient 
for computing any observable for baryon electroexcitation. Furthermore, it is suffi- 
cient to compute the matrix elements of I~^{0) in order to obtain the reduced matrix 
elements. To see this, we show that the matrix elements of the remaining components 
of P^(0) can be obtained by suitable transformations of /^(O) matrix elements. For 
spacelike momentum transfer, it is always possible to find a frame in which g"*" = 
and the spatial momentum transfer q^ lies along the x axis. Given matrix elements 
{M'j'] P'iJ,'\I~^[0)\Mj; P/i) in this frame, the matrix elements of /~(0) can be obtained 
by a rotation of tt about the x axis: 

(M'j';PV|/-(0)|Mj;P/x) = (M'/; PVf/t[P,(7r)]|J+(0)|f/[P,(7r)]Mj; P^), (28) 

and the matrix elements of /^(O) can be obtained by a rotation of the matrix element 

/3=l(J+_J-)by|: 

(MY;PV|/^(0)|Mj;P^) = (M'/; PVf/t[P,(^)]|/3(0)|f/[P,(^)]Mj; P/i), (29) 
The matrix element of /^(O) is constrained by the continuity equation: 

(M'j';PV|/i(0)|Mj;P/i) = |^(M'j'; PV|/+(0)|Mj; P/x), (30) 



where Q = \fQ^. 

4.1 The Q2 ^ limit 

Our calculation should go smoothly over to the real-photon case when we calculate in 
the limit Q^ -^ 0. Note that since q^ = 0, we have Q^ = — q^. This means that the 



spatial momentum transfer q^ vanishes as Q^ -^ 0, and as a consequence all of the 
light-front matrix elements (M'j'; P'/i'|J+(0)|Mj; P/x) vanish due to orthogonality 
between the initial and final states. In this limit, the continuity equation, Eq. (pUD, 
becomes 

(My;PV|/i(0)|Mj;P/i) = |-(My; PV|/+(0)|Mj; P/i), (31) 

and q^ = (g°,0,0,— g°) = (q'~,q± = 0,q~^ = 0), so the spatial component of q^ 
lies along the 2;-axis. This means that the matrix elements of /^(O) give the desired 
transverse-photon amplitudes. Since both the matrix elements of /^(O) and q^ are 
tending to zero, we can rewrite Eq. (RTf) 



(MY;Py|/i(0)|Mj;P/i) = lim ^ (M'j';PV|/+(0)|Mj;P/x). (32) 

q2^o 2 oq± 

In practice, rather than calculating a derivative, we use the multipole amplitudes 
to compute helicity amplitudes, and the latter have smooth behavior as Q^ ^ 0. 

Note that if current continuity is imposed on our calculation at the level of the 
light-front matrix elements, i.e., adherence to Eq. (|30D is ensured by writing the 
/^(O) matrix elements in terms of those of /"""(O) using Eq.( 0), then we will have 
a singularity at Q^ = —q]_ = 0. Implementing Eq. ( p^ ) is technically difficult when 
the /^(O) matrix elements are evaluated approximately. We have therefore chosen 
[e.g., for the A^ to A transition, which has three independent multipole amplitudes) 
to calculate three of the four multipole invariants from three of the four /"^(O) matrix 
elements, and we impose continuity at this level by writing the fourth invariant in 
terms of the other three. 

In practice we have chosen to eliminate the light-front matrix element correspond- 
ing to the largest change \fi' — fi\ in magnetic quantum number for the light-front state 
vectors in Eq. (^Tf). The choice of which multipole to constrain is then usually obvi- 
ous; for example in the case of a transition between a nucleon an excited state with 
J^ = I , there are three multipoles: the Coulomb multipole with (/, s, J7) = (0,0,0), 
magnetic (1,1,1), and electric multipoles (1,1,0). Clearly we cannot eliminate the 
magnetic multipole, since continuity relates the Coulomb and longitudinal matrix el- 
ements. Elimination of the Coulomb multipole in favor of the electric multipole leads 
to unphysical results for Ci at intermediate values of Q^ (where g*^ = 0). Therefore, 

2 

in this case we constrain the electric multipole by imposition of current continuity. 

Similarly, in the case of J-'^ = | states, we have the magnetic (/, s, J7') = (l,l,l), 
Coulomb (2,0,2), electric (1,1,2), and second electric (3,1,2) multipoles. Once again a 
physical solution cannot involve elimination of the magnetic or Coulomb multipoles; 
we have chosen to eliminate the electric multipole with highest / by imposition of 
continuity. Making the other choice changes our results for the helicity amplitudes 
by an amount comparable to or smaller than the rotational covariance condition 
described and evaluated in what follows. Similarly, for J^ = ■^ and | states we 
also eliminate the electric multipole with highest /. 

9 



4.2 Further conditions on the matrix elements 

While the matrix elements of /^(O) are sufficient to determine the reduced matrix 
elements introduced in Eq. (^2]) , they are in fact not independent of each other. Parity 
considerations imply that 

{M'f; P' - /x'|/+(0)|Mj; P - ^l) = n'n(-iy'-^(-l)^'-^(M'j'; P'/i' |/+(0)|Mj; P/i) 

(33) 
This cuts the number of independent matrix elements in half. In addition, there can be 
constraints which come from the requirement of rotational covariance of the current 
operator. To derive the constraint conditions, we note that, for matrix elements 
between states with canonical spins, in a frame where the three-momenta P' and P 
lie along the quantization {z) axis, 

,{M'f- P' - ^'|/^(0)|Mj; P - /i)e = 0, 1^' - ^1 > 1. (34) 

That is, helicity must be conserved. Since Eq. (0) must be satisfied for all compo- 
nents of the current, it must also be satisfied for matrix elements of /'^(O). Trans- 
forming to light-front momenta and spins, with momentum transfer along the x axis, 
we obtain 



5:Z}^t,,(E;j(MY;P'A'|/+(0)|Mj;PA)Di^(Ee.) = 0, |/x' - /x| > 2. (35) 



A'A 

The rotations above are 



Ech = Ec/(P, M)E,(^), BLh = EcfiP', M)Ry{^), (36) 

where R^f is the Melosh rotation of Eq. (^ which, together with the rotation Ry{^), 
transforms the state vectors from light-front spin to helicity. For elastic scattering, 
Eq. ( P5| ) is applicable only to targets with j > 1. For elastic and inelastic scattering 
involving higher spins, there is a separate, unique rotational covariance condition for 
each pair of helicities whose difference is two or more. Thus, for a transition | ^ f , 
there is a single condition, while for | — ^ f , there are three. Helicity pairs which 
differ by an overall sign change do not generate additional conditions. 

The requirement of rotational covariance provides a dynamical constraint which 
cannot be satisfied without the introduction of interaction-dependent currents, i.e., 
two-and three-body current operators. Thus, while the parity constraints in Eq. ( pBf ) 
can be satisfied in a calculation employing one-body current matrix elements, the 
rotational covariance condition in Eq. (|35|) cannot. A measure of the violation of 
the condition (and hence the need for many-body current matrix elements) is then 
the value of the left-hand side of Eq. (^) for each independent pair of helicities for 
which the condition is nontrivial. In an earlier work |jT6[ , it was shown that rotational 



covariance tends to break down for constituent models of mesons when j^ ^ 1. 
For mesons with mass of a few hundred MeV, this limits the applicability of such a 

10 



calculation to Q^ less than 1 or 2 GeV^. Since baryons are hundreds of MeV heavier 
than light mesons, we would expect the violation in this range of Q^ not to be so 
severe, but it can be checked directly from the calculated matrix elements, and is 
discussed further below. 



5 Results 



The result of combining Eqs. (|l6D -(p!9|) is a six- dimensional integral over two relative 
three momenta. These integrations are performed numerically, as the angular inte- 
grations cannot be performed analytically. The integration algorithm is the adaptive 
Monte Carlo method VEGAS [^. Typical statistical uncertainties are on the or- 
der of a few percent for the largest matrix elements. In what follows we have taken 
point-like constituent quarks, i.e., with FiqiQ"^) = F2q{Q^) = 1, in our evaluation of 
Eq.( ITBD. The light-quark mass is taken P, E3| to be m„ = m^ = 220 MeV. 



5.1 Nucleon elastic form factors 

Using the techniques outlined above we can form the light-front current matrix ele- 
ments for nucleon elastic scattering {M^ ^; P'fi'\I^{0)\MN |; P/^), from Eq.( |^). We 
have evaluated Eq.( ^) using a simple ground-state harmonic oscillator basis state, 

$o,o(A;„ K,) = ^^exp {-[k^ + Kl]/2alo]} , (37) 

where the oscillator size parameter ano is taken |21|, [T|] to be 0.41 GeV. Eq. (|T8D 
applies equally well to quark spinor and nucleon spinor current matrix elements, so 
we can extract Fi((5^) and F2{Q'^) for the nucleons directly from the above light-front 
matrix elements. 

Figure |I| compares the proton and neutron Ge and Gm calculated in this way, and 
by using the same wavefunction and the usual nonrelativistic approach. Also plotted 
in Fig. |T] is the modified-dipole fit to the data. Our (Breit-frame) nonrelativistic 
calculations use a quark mass of rriu^d = 336 MeV (from a nonrelativistic fit to the 
nucleon magnetic moments) and the same oscillator size parameter as above, with 
the first order nonrelativistic reduction of the electromagnetic interaction operator 

^" = - E {|^ (P^ ■ A. + A, ■ p,) + /i,o-, ■ B,| , (38) 

where m = rriu = rn^ is the the constituent quark mass, Cj, crj/2, and /ij = gei/2m 
are the charge, spin, and magnetic moment of the quark i, and A^ = A(rj). 

Our choice of quark mass for the relativistic calculation, while motivated by pre- 
vious work P, ^, gives a reasonable fit to the nucleon magnetic moments. The 
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Figure 1: Proton and neutron elastic form factors Ge and Gu- The solid curves 
are the relativistic calculation, the dotted lines give the corresponding nonrelativistic 
result, and the dashed lines are the modified-dipole fit to the data. 

relativistic calculation yields a proton charge radius close to that found from the 
slope near (^^=0 of the dipole fit to the data. The nonrelativistic calculation falls off 
too rapidly at larger Q^ [like exp(— q^/Ga^o)]) which is not the case for the relativistic 
calculation. These observations confirm those made earlier by Chung and Coester 0] 
and by Schlumpf 0. 

5.2 Helicity amplitudes 

For spins other than j' = | it is convenient to compare the results of our calculation 
to helicity amplitudes. These are defined in terms of the matrix elements found in 
Eq. (^21) from our multipole invariants as follows: 



A^l 



N 



C\ 



Ana 



2 



2Kw 
( 

1 
( 
2 V '^Kw 



2K^ 



w 






C^/;S^c(M^./;Pl|/ii 



Mn^; P^) 



(39) 



cJ^c(M^./;P'i|/o° 



M^l-i^i 



P|) 



2/c, 



where ( is the sign of the Ntt decay amplitude of the resonance N*, a = e'^/An 
1/137, Kw is the equivalent real-photon cm. frame three momentum Kyy = {W"^ 
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M'^)/2W, and kc.m. is the virtual-photon cm. frame three momentum 

In the above Q = \/—(p is the magnitude of the four-momentum transfer, and W 
is the square root of the invariant mass evaluated at resonance, where W = Mpf*. 
Note that the square root factors are introduced in order that the Q^ —>■ limit of 
the electroproduction amplitudes corresponds to the photocoupling amplitude; the 
only restriction on these factors is that their limit when Q^ = is the same as 
the normalization factor for an external real photon. The above choices, therefore, 
represent a convention. Other authors calculate the quantity 5*1 = |kc.m. |C'i/Q- Note 

2 2 

we use an operator /°(0) for the 0-component of the electromagnetic current which is 
defined so that c(M^|; P'||/o(0)|M7v|; P|)c = +1 at Q^ = q, where iV is a proton. 

5.2.1 Helicity amplitudes for A(1232) electroexcitation 

Figure || shows our results for the Ai, A3, and Ci helicity amplitudes for elec- 

22 2 

troexcitation of the A| (1232) from nucleon targets, compared to the nonrelativistic 
results. In this case we have used the simple single oscillator-basis state from Eq. (|37D 
for both the initial and final momentum-space wavef unctions. The parameters ano 
and rUu^d are the same as above. Although the relativistic calculation does not solve 
the problem of the long-standing discrepancy between the measured and predicted 
photocouplings, the behavior of the relativistic calculation is closer to the faster- 
than-dipole fall off found in the data. The data show no evidence for the initial rapid 
rise with Q^ shown by the nonrelativistic calculation, as pointed out by Foster and 
Hughes 13. 

We have also plotted the numerical value of the rotational covariance condition 



(multiplied by the normalization factor (J47ra/2Kw for ease of comparison to the 

physical amplitudes), given by the left-hand side of Eq. (^5]), for |/i' — /i| = 2. At 

lower values of Q^ the rotational covariance condition expectation value is a small 

fraction of the transverse helicity amplitudes, but approximately the same size as Ci 

2 
and larger than the value of E2/M1 implied by our Ai and A3. Calculations which 

2 2 

attempt to predict the ratio using this approach |18[ will in general be limited by 



rotational covariance uncertainties of similar magnitude. 

5.2.2 Helicity amplitudes for electroexcitation of radially excited states 

Given the controversy surrounding the nature of the baryon states assigned to radial 



excitations of the nucleon and A(1232) in the nonrelativistic model |2^, we com- 
pare nonrelativistic and relativistic calculations for simple basis states which can be 
used to represent the two Pn resonances, A^(1440)| and A^(1710)| , as well as the 
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Figure 2: A(1232) electroexcitation helicity amplitudes A^% (solid curve) and A^% 
(dashed curve), and Cf'J^ (long dashes). The dotted curves give the corresponding 
nonrelativistic results for the Breit-frame transverse helicity amplitudes {Cf'J^ = 
in the nonrelativistic model). Also plotted is the value of the rotational covariance 
condition (dot-dashed curve), which happens to nearly coincide with Cy^. 



A(1600)| , which is assigned as a radial excitation of A(1232) in the nonrelativistic 
model. 

Figure ^ shows our results for the Roper resonance, N^ (1440), for both proton 
and neutron targets. Here we have used the single oscillator-basis state Eq. (|37D for 
the initial momentum-space wavefunctions, and a radially excited basis state 



(!>^'n(k,,Kx) 



{K 



■0,OV"'P? 



\/2,Ti2aljo 



?>-[kl + Kl]/alo) exp {-[^J + Kl]/2alo\} (41) 



is used to represent the excited final state. The sign C, of the Nn decays amplitude used 
here is calculated in the ^Pq model of Ref . , using exactly the same wavefunctions. 
As a consequence, the sign of our nonrelativistic Roper resonance photocouplings 
calculation [|l^] disagrees with the sign which appears in Ref. 0, where a reduced 
matrix element was fit to the sign of the Roper resonance photocouplings. 

There are large relativistic effects, with differences between the relativistic and 
nonrelativistic calculations of factors of three or four. Interestingly, the transverse 
amplitudes also change sign at low Q^ values approaching the photon point. The 
large amplitudes at moderate Q^ predicted by the nonrelativistic model (which are 
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Figure 3: A^(1440) electroexcitation helicity amplitudes A^^ and A"^ (dashed 
curves), and Cf,2 and C]y2 (long-dashed curves). The dotted curves give the cor- 
responding nonrelativistic result for the Breit-frame transverse helicity amplitudes, 
and the dot-dashed curve gives the nonrelativistic C^i^ [Cyi^ is zero in the nonrela- 
tivistic model). 



disfavored by analyses of the available single-pion electroproduction data [^) appear 
to be an artifact of the nonrelativistic approximation. This disagreement, and that 
of the nonrelativistic photocouplings with those extracted from the data for this 
state [O] , have been taken as evidence that the Roper resonance may not be a simple 



radial excitation of the quark degrees of freedom but may contain excited glue |2^, 
^ . The strong sensitivity to relativistic effects demonstrated here suggests that this 
discrepancy for the Roper resonance amplitudes has a number of possible sources, 
including relativistic effects. 

We also find in the case of proton targets that there is a sizeable C\ , reaching a 

2 
1 

value of about 40 x 10^'^ GeVs at Q^ values between 0.25 and 0.50 GeV^. Correspond- 
ingly, there will be a sizeable longitudinal excitation amplitude. The nonrelativistic 

C^ amplitudes shown here and in all other figures are listed (up to the A^vr sign Q 

2 
in Appendix C; formulae for the nonrelativistic transverse amplitudes (up to the sign 

C) are tabulated in Ref. |1|]. 

In the case of the N{\1VS)\ , which is assigned a second radially excited wave- 
function in the nonrelativistic model, our results using a simple wavefunction made 
up of a linear combination of basis functions (see Refs. ||21[ and [0]) are shown in 
Fig. ^. Here we see an even greater reduction in the size of the transverse ampli- 
tudes in the relativistic calculation. The C1/2 amplitudes at the photon point also 
are greatly reduced in size, although at larger Q^ values they become comparable to 
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Figure 4: A^(1710 electroexcitation helicity amplitudes. Key as in Fig. 3, except here 
the nonrelativistic 0^2 (dot-dashed curve) is nonzero. 

the nonrelativistic charge helicity amplitudes, with a similar modification of the Q^ 
dependence to that shown above in the case of the Roper resonance and G^. 

It is interesting to see whether this pattern is maintained in the case of the 
A(1600)| , which has a similar relationship to A(1232) as the Roper resonance 
has to the nucleon. This is a quark-spin | state [with spatial wavefunction given 
by Eq. (|4lD] in the nonrelativistic model, and so fundamentally differs from the Pn 
states described above. However, examination of the nonrelativistic calculation of the 
electroexcitation amplitudes shows that in both cases it is the spin-fiip part of the 
0{p/m) electromagnetic Hamiltonian which is responsible for the transverse transi- 
tion amplitudes. 

This similarity persists in our relativistic calculation (Figure ^, where we see 
photocouplings which have changed sign in comparison to the nonrelativistic results, 
along with substantially reduced transverse amplitudes at intermediate Q^ values. 
The falloff with Q^ at higher Q^ values is quite gradual. As in the case of the 
A(1232), the relativistic C1/2 amplitudes (zero in the nonrelativistic model) are small 
and comparable in size to the rotational covariance condition, which is generally a 
small fraction of the transverse amplitudes. 

5.2.3 Helicity amplitudes for electroexcitation of P-wave baryons 

We have also calculated helicity amplitudes for the final states A^| (1520) and N^ (1535), 
for both proton and neutron targets. Here we use the same initial-state wavefunction 
as above and final state wavefunctions which are made up from the orbitally-excited 
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Figure 5: A(1600) electroexcitation helicity amplitudes. Key as in Fig. 2 (here also 
the nonrelativistic amplitudes C^?2 are zero). 
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(42) 



(where /c_,_ = fc^. + iky) and their counterparts reached by angular momentum lower- 
ing. In this case configuration mixing due to the hyperfine interaction is included in 
the final-state wavefunctions. Since these states are degenerate in mass before the 
application of spin-spin interactions, they are substantially mixed by them; details 
are given in Appendix B. 

The results for the helicity amplitudes for A^| (1535) excitation from both pro- 
ton and neutron targets are compared to the corresponding nonrelativistic results in 
Figure ||. 

In contrast to the results shown above, in this case there appears to be little 
sensitivity to relativistic effects in the results for the transverse amplitudes A1/2; this 
is not the case for the C1/2 amplitudes. For both targets there are substantial C1/2 
amplitudes at small Q^, and the Q^ dependence is very different in the relativistic 
calculation (resembling the dipole behavior of the nucleon form factors in Fig. ||). 
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Figure 6: A^(1535) electroexcitation helicity amplitudes A^/g and y4"/2 (dashed 
curves), and Cj'/g and Cjyg (long-dashed curves). The dotted and dot-dashed curves 
give the corresponding nonrelativistic result for the Breit-frame transverse and C1/2 
helicity amplitudes, respectively. 

Our results for the A\ and C\ helicity amplitudes for A^| (1520) excitation from 
2 2 

both proton and neutron targets are compared to the nonrelativistic results in Fig- 
ure ^ and this comparison for Az and the value of the rotational covariance condition 

2 
are shown in Figure ^. 

Here we see large relativistic effects in the A\ amplitudes. From Ref. |[l| one can 

2 
see that the A\ amplitudes for this state are proportional to 1 — k"^ /o^jq for a proton 

2 
target, and 1 — k"^ / {"ialjo) fo^ ^ neutron target, where k = |k| is the virtual-photon 
three-momentum. Here the constant term arises from the convection part of Hem, and 
the k'^ /oi^HO term arises from the quark-spin- flip part. Our relativistic treatment can 
be expected to change the relationship between these two terms (as well as adding 
other effects), and we see here that this has caused a substantial cancellation. In 

this case the charge amplitudes Ci are greater in magnitude than the transverse 

2 
amplitudes of both proton and neutron targets, and they are the largest at small Q^ . 

From Fig. ^we can see that the relativistic effects in Az are smaller (in the nonrel- 

2 
ativistic model there are not two partially cancelling terms), except near Q^ = 0. In 

the case of a proton target, the rotational covariance condition expectation value is a 

substantial fraction of the A 3 amplitude at larger values of Q^, and is comparable to 
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Figure 7: A^(1520) electroexcitation helicity-^ amplitudes; key as in Fig. 6. 
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the small A^i at all Q^ values; the absolute size of the rotational covariance condition 

2 

can be taken to be a measure of the absolute uncertainty introduced in our results 
from considerations of rotational covariance. 

As an example of a predominantly quark-spin- 1 P-wave excitation, Fig. § shows 
relativistic and nonrelativistic amplitudes for electroexcitation of N{1650)^ , the 
partner of the predominantly quark-spin-^ state A^(1535)| (see Eqs. ([ 
and (^) in Appendix B). 
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Figure 9: A^(1650) electroexcitation helicity amplitudes. Key as in Fig. 6. 

Here, in contrast to those of iV(1535)| , the transverse helicity amplitudes show 

considerable sensitivity to relativistic effects near Q^ = for proton targets, and 

at all Q^ values for neutron targets, and as above the charge helicity amplitudes 

have remarkably different Q^ behavior. Amplitudes for the predominantly quark- 

spin-| state A^(1700)| tend to be quite small if calculated in the nonrelativistic 

model [^ |12[, and this is still true in our relativistic calculation, with the exception 

of the amplitude A^ . 
2 
We complete our survey of relativistic effects for basis functions representative of 

P-wave states by examining the amplitudes for electroexcitation of A(1620)| , once 

again using a simple linear single-oscillator basis wavef unction [Eq. ( |53| ) in Appendix 

B]. Fig. contrasts our results with the nonrelativistic results of Ref. [|21[ for Ai and 

2 
for Ci from Appendix C In both cases the amplitudes near Q^ = Q are considerably 
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Figure 10: A(1620) electroexcitation helicity amplitudes. Key as in Fig. 6. 

larger in the relativistic calculation, with the transverse amplitude decreasing to below 
the nonrelativistic result above approximately Q^ = 0.7 GeV^ due to quite different 
(dipole-like) behavior as a function of Q^. Note in the case of negative-parity A 
states, the charge helicity amplitudes are not zero if calculated nonrelativistically 
(see Appendix C). We have found similar relativistic sensitivity in a calculation using 
a simple basis function which can be used to represent the state A(1700)| . 

6 Discussion and Summary 

The results outlined above establish that there will be considerable relativistic effects 
at all values of Q^ in the electroexcitation amplitudes of baryon resonances, even at 
Q^ = 0. In the real-photon case, this can be understood in terms of the sizeable 
photon energy required to photoproduce these resonances from the nucleon, which 
implies a photon three-momentum comparable to the quark Fermi momentum in the 
nucleon. In particular, our results show that the Q^ dependence of the nonrelativistic 
amplitudes is generally modified into one resembling a dipole falloff behavior, as has 
been shown in the case of the nucleon form factors. This behavior is likely to be 
partially modified by the inclusion in the wavefunctions of mixings to higher shells, 
which are required by any model of baryon structure which takes into account the 
anharmonic nature of the confining potential between quarks. However, we consider 
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it remarkable that relativistic effects account for a large part of discrepancy between 
the nonrelativistic model's predictions and the physical situation. 

We wish to stress that the comparisons we have made here, using simple basis 
functions representing the individual states, are intended to be representative of the 
degree of sensitivity to relativistic effects for states of their quantum numbers. In 
order to make reliable predictions for these amplitudes we should use solutions of the 
relativistic three-body problem; at the very least, configuration mixing of the kind 
present in the nonrelativistic model must be included in both the nucleon and all 
of the final states, as it has been shown to have substantial effects on the predicted 
amplitudes |1|, |10|, |12|. It is for this reason that we have not compared our results 



directly to the limited available data 0. 

We also note that the Hamiltonian used in Ref. |^ is in fact a three-body mass 
operator, and so its eigenfunctions can be used directly and consistently in a calcu- 
lation of the kind described in this paper. We have developed better techniques for 
performing the numerical integrals required to obtain the light-front matrix elements 
between configuration-mixed states, whose more complicated integrand structures 
make standard Monte Carlo methods inefficient. 

Nevertheless, it is obvious from the results presented above for radially excited 
basis states that electroexcitation amplitudes of the Pn Roper resonance A^(1440)| 
and A^(1710)| states, as well as those of the A(1600)| , are substantially modified 
in a relativistic calculation. Given the controversial nature of these states 1^^ |28|, we 
consider this an important result. Our results show that relativistic effects tend to 
reduce the predicted size of the amplitudes for such states at intermediate and high 
Q^ values, in keeping with the limited experimental observations for the best known 
of these states, iV(1440)i^. 

Within the context of a model such as this one, full gauge invariance cannot be 
achieved. Reasonable results can be obtained by constraining one of the higher multi- 
poles in terms of lower multipoles. This can be thought of as a variant of the Siegert 



hypothesis pOf, which has been used successfully in the study of electromagnetic 
properties of nuclei for many years. 

We have also found that the rotational covariance violation is a small fraction of 
the larger amplitudes for the Q^ values considered here. In cases where the dynamics 
causes an amplitude to be intrinsically small, the uncertainty in our results for these 
amplitudes becomes larger. In particular, the calculated ratios E2/M1 and C2/M1 
for the electroexcitation of the A(1232)| in the absence of configuration mixing of 



D-wave components into the initial and final state wavefunctions |T3] are probably 
100% uncertain, and are thus consistent with zero at all Q^ |^ . This may not be the 
case in the presence of such configuration mixing, and we intend to investigate this 
possibility, since A(1232) electroproduction is the subject of current experiments at 



MIT/Bates and several proposed experiments at CEBAF |^ 



22 



Acknowledgements 

We would like to thank Mr. Heng Rao for his help with the analytic nonrelativis- 
tic calculations. We also thank Mr. Mark Schroeder for his help with developing 
message-passing code for parallel architecture machines. The calculations reported 
on here were made possible by grants of resources from the Pittsburgh Supercom- 
puting Center and the National Center for Supercomputing Applications. The work 
was also partially supported by the National Science Foundation under grant PHY- 
9023586 (B.D.K), by the U.S. Department of Energy through Contract No. DE- 
AC05-84ER40150 and Contract No. DE-FG05-86ER40273, and by the Florida State 
University Supercomputer Computations Research Institute which is partially funded 
by the Department of Energy through Contract No. DE-FC05-85ER250000 (S.C). 

Appendix A: Conventions for Baryon 
Electroproduction amplitudes 

In the following we have collected the information required to translate between stan- 
dard conventions for (pion) electroproduction amphtudes and our calculated helicity 
amplitudes. 

A.l Pion electroproduction helicity elements 

In the following we relate the helicity elements Ai±, Bi±, and Ci± conventionally 
used to describe pion electroproduction to the helicity amplitudes defined above in 
Eq. (pQl). The Ai±{Bi±) elements are present when a resonance is excited by photons 
with their spin projections antiparallel (parallel) to the target-nucleon spin projection; 
the momentum of the incoming photon is parallel to the direction of the initial nucleon 
spin. The Ci± elements are present when a longitudinal photon excites the resonance 
(there is no change in the target nucleon spin projection). These can be written in 
terms of the helicity amplitudes as 

2 



^l± — T^TrNJ^l ) 
2 

where C^^ is an isospin Clebsch-Gordan coefficient for A^* -^ Ntc, J is the total an- 
gular momentum of the resonance A^*, and |fcc.m.| is the cm. frame three-momentum 
of the photon. The / in the subscripts above is the angular momentum of the relative 
wavefunction of the final-state A^ and vr. Note that the J of the resonance must be 
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J = / ± i, so for a given J there are two possible / values. However only one of 
these values is consistent with parity conservation. For example for the A(1232) with 
J^ = I we have / = 1 or / = 2, but a nucleon (positive intrinsic parity) and a 
pion (negative intrinsic parity) in an / = 2 relative wavefunction have negative overall 
parity, so this value is not allowed. The sign in the subscripts above is positive if 
J = Z + i and negative ii J = I — \. So for the A(1232) only Ai+, i?i+ and Ci+ are 
possible. 

The Breit-Wigner factor / for the decay A^* 



A^TT at resonance is 



/ 



K^ 



w 



MmT 



N ^ Nn 



(2J+l)7r|(g; 



W 1.2 



tot 



where (g7r)c.m. is the pion cm. frame three-momentum 

[W^ _ (^^ + M^fW - (m^ - M^)2] 



7? 
i''"/c.m. 



41^2 



(44) 



(45) 



With these conventions the A^7 partial width for a resonance to decay to nucleon N 
is 
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(46) 



A. 2 Electric, magnetic and scalar multipoles 

In the following we relate the magnetic, electric, and scalar/longitudinal multipoles 
Mi±, Ei±, and Si± to the hehcity elements described above; these can then be related 
to the helicity amphtudes calculated here by use of Eq. (^3|). A transition is defined to 
be magnetic with multipoles Mi± if the total angular momentum absorbed from the 
photon is /, and electric with multipoles Ei± if the total angular momentum absorbed 
from the photon is / ± 1. The corresponding scalar /longitudinal multipole is Si±. If 
J = Z + I the scattering can be described in terms of the magnetic, electric and scalar 
multipole amplitudes Mi^, Ei^, and 5"/+, which can be written in terms of Ai^, Bi^ 
and C/+ as 



Ei+ 



1 



2(/ + 
1 



- [2Ai+ -{1 + 2] 



B, 



i+\ 



S, 



1+ 



2(/ + l 
1 \k 

J^ I "'cm 

/ + ! Q 



[2A 



1+ 



IB, 



i+\ 



(47) 



^C; 



/+• 



If J = / — ^ the scattering is described by Mi^, Ei^, and Si^ with 



Mi_ = -[2Ai_ + {l-l)Bi_ 



24 



El. 



Si. 



- [-2A- + (/ + l)5z-] 



(4J 



1 \kr 



The inverse transformation is 



/ Q 



-C, 



Ai+ = -[lMi+ + (l + 2)Ei+] 

Bi+ = -M1+ + E1+ 

A^- = i [(/ + 1)M,_ - (/ - 1)E,_ 

Bi_ = M_+Ei_. 



(49) 



A. 3 Invariant multipoles for A(1232) electroproduction 



For the special case of electroproduction of a J^ 



particle such as the A I (1232), 



we can define invariant multipoles G%j, G% and G^. The relation between these and 
the helicity amplitudes defined above is 



E 

G*c 



-F 



-F 



V3A^ + A^ 
2 2, 

^Af - A^ 
v3 2 2 



2^FV2Sf, 

rhn m ^ 



(50) 



where 



Note that Eq. 
Kw) we have 



1 M 



N 



ATia W 



K 



w 



1 + 



Q' 



{Mn + Wf 



2 M, 



Af 



(51) 



]) implies that at the photoproduction point {Q^ = 0, where \kc 



F(0) 



Ana\ 



2M% 



W^-Ml' 



(52) 



Appendix B: Wavefunctions for the P-wave 
baryons 

We construct A^ and A states which have the correct permutational symmetry (mixed- 
A type for A^ states, and symmetric for A states) in their flavor/spin/space wave- 
functions. A common totally antisymmetric color wavefunction is assumed, as is an 
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implicit Clebsch-Gordan sum coupling L and S. The spin-quartet A^ states (notation 
is \F'^^^^L(jJ^), with F the flavor A^ or A, and a the spatial permutational symme- 
try) \N^Pm\ ), \N^Pm^ ), and \N^Pm^ ) are made up of Clebsch-Gordan linear 
combinations of 

uudxi ^iM ■ (53) 

The spin-doublet A^ states \N'^Pm\ ) and \N^Pm^ ) are made from 

uud^(x1 Km-x\ ^IaX (54) 

and the spin-doublet A states \A'^Pm^ ) and |A^Pjv/| ) are made from 

uud^lx'i Km + x\ ^i.mV (55) 

Where mixing is allowed, the physical states do not have definite quark spin S. We 



have used 21 



|iV(1535)|") ^ 0.85\N^PMf) + 0.53\N''PMf) 

\N{1650)1') -^ -0.53\N^PMf) + 0.85\N^PMl') 

|Ar(1520)|") ^ 0.99\N^PMf) ~0.n\N^PMf) (56) 

|iV(1700)|") ^ 0.11\N^PmI') + 0.99\N^PmI'). (57) 

Appendix C: Nonrelativistic Ci amplitudes for 

2 

electroexcitation of resonances 

In the Figures above we have plotted predictions of the Isgur-Karl-Koniuk (nonrel- 
ativistic) model for the charge helicity amplitudes C\ and C". To the best of our 

2 2 

knowledge these are not published elsewhere, so we list them below. 
For Ari^(1440), 



2 



^/Sk'^niq 



9a^ 



HO 



:i,o), (5^ 



where 



Ao = #;T^e-^'^/(^"-), (59) 

V fco 2mg 

and (A;o,k), with k = |k|, is the (virtual) photon four-momentum in the frame in 
which the calculation is performed (taken here to be the Breit frame). For the states 
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and N^Pm^^ 


\ J^XJ'Ljyj J n VVV^^1VV_/ UAiV_/ 


dilJ-Aj^/llU C 


iVJ-V^O L\J1. U 


LiV^lX U VVV^ "^ 




2 


= Ao 


.V^krriq 
3aHo 


(1,-1), 




2 


= 0, 


and similarly for the states A^(1520)|~ 


andiV(1700)|" 




c^r\N'PMf) 

2 


= Ao 


2kmg 
i 

Sano 


(-1,1), 




c^r\N'PMf) 

2 


= 0. 


For Ai"(1620) 











PI- 



and for Af (1700) 



For iVl^(1710), 



Cl^ = -Aot 



V2k 



m„ 



iN 



Sano 
2km„ 



2 oauo 



C'l 



{p,n) 



An 



V6k' 



m„ 



l^ajjo 



:i,-i). 



(60) 



(61) 



(62) 



(63) 



(64) 
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